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Abstract 

Two conditions are derived for Ising models to show non-universal critical 
behaviour, namely conditions concerning 1) logarithmic singularity of the 
specific heat and 2) degeneracy of the ground state. These conditions are 
satisfied with the eight-vertex model, the Ashkin- Teller model, some Ising 
models with short- or long-range interactions and even Ising systems without 
the translational or the rotational invariance. 
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§ 1 Introduction 



The universality of critical exponents is one of the most important concepts 
in critical phenomena. According to this universality hypothesis, critical ex- 
ponents depend only upon the dimensionality, the symmetry and the inter- 
action range of Hamiltonians, namely they are independent on the details of 
the relevant Hamiltonian such as the strength of the interactions in ordinary 
situations. 

It is quite interesting to study in what condition this hypothesis is violated 
from the form of the Hamiltonian. 

The first example violating the universality hypothesis is the eight-vertex 
model solved by Baxter [1]. This model can be mapped[2][3] onto the two- 
layered square-lattice Ising model in which two layers interact each other via 
a four-body interaction J4. The critical exponents of this model vary with the 
parameter \x which is a function of interaction energies. The exponents are 
obtained[l] as a = 2 — n/fi, f3 = n/16fi, v = n/2fx, and the scaling hypothesis 
or the weak universality hypothesis[4] insists that 7 = 7u/A = 7^/8^, 5 = 15, 
77 = 1/4. Kadanoff and Wegner[3] have shown that the existence of marginal 
operators is a necessary condition for appearance of continuously varying 
critical exponents. Kadanoff and Brown [5] have shown that the long-range 
behaviour of the correlations of the eight- vertex and the Ashkin- Teller models 
are asymptotically the same as those of the Gaussian model. 

There exists another model which consists of short-range two-body inter- 
actions and which is believed to have continuously varying critical exponents. 
The s = 1/2 square-lattice Ising model with the nearest-neighbour interac- 
tion J and the next-nearest neighbour interaction J' have been studied by 
many authors [6] -[12] to obtain the phase diagram. The numerical calcula- 
tions of van Leeuwen,[13] Nightingale [14] and Swendsen and Krinsky[15] are 
the first to show non-universal critical behaviour of this model. The singu- 
lar part of the free energy is calculated perturbatively by Barber [16]. The 
results of the high temperature expansion by Oitmaa[17] agrees with these 
calculations. The coherent-anomaly method (CAM) is applied[18]-[22] to this 
model and the continuously varying critical exponents are estimated with er- 
rors smaller than ~ 1%[21][22]. The degeneracy of the ground state energy 
and the existence of a multi-component order parameter have been studied 
by Jiingling[23] and by Krinsky and Muhamel[24]. A more general Hamil- 
tonian has been investigated by the present authors [22] which includes the 
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above two models as special cases of it and which has continuously varying 
critical exponents. 

In the present paper, we phenomenologically derive a sufficient condition 
to have continuously varying critical exponents. This study is the general- 
ization of the argument reported in ref. [25]. Our condition is satisfied in the 
eight- vertex model, the Ashkin- Teller model and the s = 1/2 square-lattice 
Ising model with the next-nearest-neighbour interaction. Our condition is 
also satisfied with some systems including long-range interactions and with 
some systems without the translational or the rotational invariance. 

Brief explanations of the relevant models are given in section 2, together 
with some explicit conditions on continuously varying critical exponents. A 
phenomenological perturbation scheme [26] is explained in section 3 and ap- 
plied to the eight-vertex model in order to demonstrate that our scheme 
provides the exact first-order derivative of the critical temperature and the 
exponent 7 of the eight-vertex model in section 4. The temperature depen- 
dence of correlation functions and the characteristic cancellations of interac- 
tion energies at the ground state are discussed in section 5. The symmetry 
of the relevant model is studied in section 6. The perturbational scheme 
explained in section 3 is applied, in section 7, to the model defined in section 
2. It is derived that there exist finite first- or second-order derivatives (with 
respect to the interaction energies) of the critical exponent 7 and hence it 
varies continuouly as a function of interaction energies. Finally, the condition 
in section 2 is generalized to include a more general type of interactions in 
section 8. 



§ 2 Models 

Let us consider the following Hamiltonian 

k kl 

where {Hk} is a finite set of two-dimensional Ising systems. In the present 
paper, we derive that the model with the Hamiltonian (2.1) have continuously 
varying critical exponents provided it satisfies the following two conditions 
1) and 2). 
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1) The specific heat of Tik shows the logarithmic singularity at the critical 

temperature T c which is independent on k. 

2) The ground state energy of the total Hamiltonian 7i is invariant for the 

spin inversion of each Ti^. Here Tiki is written as 

Wh = -JhE°W ) » ( 2 - 2 ) 

i 

where Of^ and Of* are the n( fc )-body and n^-body spin product of 
spins belonging to TCk and Hi, respectively. 

The logarithmic singularity in the condition 1) results in the temperature 
dependence of the correlation functions of the form (siSj) = Cq + Cielne, 
where e = (T — T c )/T c and Co and c\ are some constants. This formula is 
derived in section 5. 

The restriction on and Of 1 in the condition 2) can be partially 
removed and more complicated interactions are permitted, i.e. and 
can depend on the region lZ m in (5.6). This generalization is explained in 
section 8. 

§ 3 Perturbational scheme 

To derive non-universal critical behaviour of these models, let us consider 
the perturbational expansion [26] with respect to interaction energies. The 
susceptibility x is assumed to behave as 

X - e(J)-^ J \ e(J) = (T — T C (J))/T C (J), (3.1) 

where J is an interaction energy. Differentiating it with respect to J, we 
obtain 

(| )j ^ W - T( 0)^(|l) -(|l) lo ge (0)], (3.2) 

where the subscript denotes J = 0. Hence we can obtain (d~f/dJ)o from 
the coefficient of the temperature dependence Xo l°g e - On the other hand, 
the susceptibility is expressed by the two-spin correlation functions in the 
form 

X = /5a*bS^ojo< s <o s jo)> ( 3 - 3 ) 

iojo 
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where (3 = 1/ ksT, /zb is the Bohr magneton and {gi j } denote the signs com- 
ing from the emerging order. We differentiate (3.3) and estimate (d'y/dJ)® by 
comparing with (3.2). The existence of a finite and non- vanishing derivative 
of the exponent 7 is an evidence for continuously varying critical exponents 
to appear. 



§ 4 Example 

Here we give two important examples. Let us consider the zero-field eight- 
vertex model. The Hamiltonian H$v of this model is written as Hsv = 
Hi + H 2 + where Hi and H 2 denote the following Hamiltonians 

Hi = —J X (SijSi+lj+l + SijSi+lj-l), 
i+j=even 

H 2 = -J' (sijSi+ij+i + SijSi+ij-i), (4.1) 

i+j=odd 

and H12 is 

H12 = —J& Sj-jSj+ij+iSi+i-jSj-j+i. (4-2) 

This model decouples into the two square-lattice Ising models Hi and H 2 
when J4 = 0. The interaction H12 has even symmetry for the spin inversion 
of each subsystem Hi and H 2 and hence the energy of Hsv is four-fold 
degenerated in the whole temperature region. This model obviously satisfies 
the conditions 1) and 2) in section 2. The weight g io j equal 1 for J' > 0. 
Differentiating (3.3) with respect to J 4 and using the fact that the model 
decouples into two layers when J 4 = 0, we obtain the following expression 

s i jo s ij s kl s k+ll+l)o — (Si j Sij)o{ s kl s k+ll+l)o\{ s kl+l s k+ll)o, 

4 iojoijkl 

(4.3) 

where ( ) denotes the expectation value for J 4 = 0. This expression coincides 
with 
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Here ooo denotes the nearest-neighbour two-spin correlation function which 
shows the following behaviour 

1 4 7' 

^Vf + S^f (0)log£(0) (45) 

at the critical point. Comparing (4.4) and (4.5) with (3.2), we obtain 

dT T c (0) d 1 47(0) fAR . 

W J4 =° = vu> and W J4 =° = -^fwr (4 - 6) 

which are identical with the derivatives obtained from the exact result by 
Baxter. 

The square-lattice Ising model with antiferromagnetic next-nearest-neighbour 
interactions is another example. The Hamiltonian of this model is given in 
the form Hsaf = Hi + H 2 + H.12, where H.\ and Hi are given by (4.1) for 
J' < and H12 is 

H 12 = -Jj^isijSi+ij + SijSij+i). (4.7) 

The ground state of this model is ordered as the Neel state in each sublattice 
for the interaction region \J/J'\ < 2, in which the ground state energy is 
invariant for the spin inversion of each sublattice. It is in this interaction 
region \J/J'\ < 2 that this model is considered to have continuously varying 
critical exponents. This model also satisfies the condition in section 2. It is 
derived in ref. [16] and [25] that this model has continuously varying critical 
exponents with the derivatives (d^/dJ)o = and (<9 2 7 /<9J 2 )o 7^ 0. 



§ 5 Preliminary formulas 

The temperature dependence of correlation functions is specified from the 
condition 1). For the purpose to obtain multi-spin correlation functions, let 
us consider Tik = Hk — JO, where O denotes a certain spin product. The free 
energy f ks of this Hamiltonian is differentiate with respect to J. Assuming 
that f ks shows logarithmic or power-law behaviour for J = 0, it should have 
the form 

1 — f ~ a2 

fks = C ie 2 - ai loge + C 2 e 2 , (5.1) 

a 2 
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where C\ and C2 are some constants, ot\ and «2 are functions of J with 
aii(J = 0) = and a^J) — > for J — > 0. The second term converges to 
C 2 e 2 loge when J — > 0. From (5.1), 

(O) = \{%)j=o oc26(^) lo g 6 + O(6), (5.2) 

where contributions vanishing faster than e log e (for e — > 0) are included in 
0(e). The regular part of the free energy yields a certain constant term. Then 
we generally obtain the temperature dependence of correlation functions in 
the form 

(O) = c + Cl e\oge. (5.3) 

The condition 2) results in the strong cancellations of interaction energies 
at the ground state. Here we introduce notations for the ground-state con- 
figuration of the system as {g,j} and the spin product Oi = ■ ■ ■ Sj n in Tiki 
at the ground state as Gi — • • • 9i n - The ^/-dependent part of the ground 
state energy, €g(Jm), is written as 

ec(4) = -4Eft- ( 5 - 4 ) 

i 

In the case where J2i Oi has even symmetry for the spin invergion of Tik 
and 7ii (i-e. both and are even), the condition 2) is automatically 
satisfied. Otherwise, the condition 2) yields 

-Jh E^ = - J *< E(-&) thatis E ft = °- ( 5 - 5 ) 

i i i 

We exclude, from our arguments, the case that the condition 2) is asymp- 
totically satisfied only in the thermodynamic limit. Then we can rewrite the 
condition (5.5) as 

Eft = EEft- and Eft = °> ( 5 - 6 ) 

where {TZ m } denote a set of finite regions containing a finite number of spins 
and the symbol i G TZ m denotes that all the spins in Qi are included in TZ m . 
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§ 6 Symmetries and the vanishing derivatives 



6.1 Basic symmetries 

In this section, we derive some properties which depend only on the symmetry 
of the relevant model. We consider here the case N = 2 i.e., the following 
Hamiltonian Tt = Hi + H 2 + H12 for e > 0. 

We assume that Tin has odd symmetry for the spin invergion of Ti 2 . 
Then the interaction 7i 12 has the form 

<H l2 = -JY.O?Of\ (6.1) 

(2) 

where 0\ changes its sign for the spin inversion of 7^2- Hereafter we omit 
the subscript of J\ 2 and write J\ 2 as J for simplicity. 

Let us consider the transformation of spin configurations in which all the 
spins on Ti 2 are reversed and the spins on Tii are unchanged. The energy of 
the model 7i with the interaction J for the spin configuration C equals to 
the energy of the model H. with the interaction — J for the transformed spin 
configuration C . As a result, we can find one-to-one correspondence of the 
Boltzmann factor coming from the configuration C and C. Summing over 
all configurations, we find that the partition functions Z(J) and Z{— J) are 
the same 

Z(J) = Z{-J). (6.2) 

Since we consider the partition function without an external field, we obtain 
from (6.2) for the critical temperature T C (J) and the exponent a(J) of the 
specific heat as 

T C (J) = T c (-J) and a(J) = a(-J). (6.3) 

From (6.3) we can conclude that any odd derivatives of T C (J) and a(J) are 
equal to zero when they exist. 

Next we study the symmetry of the correlation functions 0Oi o j o (J) which 
are the expectation values of the following two-spin product with the weight 
corresponding the emerging order as 

- ^^g^tm „ e M -re(J)^}, (6.4) 
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and also we study the symmetry of the susceptibility 

x(J) = T,"io jo (J)~<J)- j{J \ (6-5) 

iojo 

where {gi j } are the sign corresponding to the emerging order and r is the 
distance between the site io and j . We can assume without loss of generality 
that the site io lies on 7Yi- Corresponding to the change of the interaction 
J to — J, the quantity g io j changes its sign when the site j belongs to TL2 
and then gi j Si Sj is even for the spin invergion of 7Y 2 . As a result, we again 
arrive at 

^iojo( J ) =^ojo(- J ) and X(J)=X(~J), ( 6 - 6 ) 

and hence 

u(J) = iy(-J) and 7 (J) = 7(-J). (6.7) 



6.2 Vanishing cases 

The above argument cannot exclude the case that the critical temperature 
and exponents do not differentiable at J = 0. The critical coefficient, indeed, 
behaves as a cusp and do not differentiable at J = when the system is not 
translationally invariant. Here we show using the scaling hypothesis that the 
first-order derivatives of the critical temperature and exponents exist and 
they are vanishing. Let us consider the free energy / without the external 
field as 

/~Ce 2 - a loge, (6.8) 

where C is some constant. Differentiating (6.8) with respect to the interaction 
J and taking a limit J — > 0, we obtain 



-) oe \o g e + Ce -(- 



and as a result 



1 (tf f)n 

+ C e 2 log6[(2 -«)-(— ) -(_) l O ge] (6.9) 



, dC . .da. , . de . . 

— ) = 0, (— ) = and (— ) = 0. (6.10) 
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Next let us consider the free energy / with the external field h. We 
assume the scaling form of / as 



fj(X"€,X"h)=Xfj(e,h), (6.11) 

where A is some parameter and p and q are numbers independent on e and 
h. They are related to the exponents a and 7 as 

1 , 2o - 1 . . 

a = 2 - - and 7 = — . 6.12 

p p 

From (6.10) and (6.12), we obtain (dp/dJ)o = 0. The partial derivative of / 
shows the same scaling form as 

dfj(\>e,\ q h) = x dfj(e,h) 
dJ dJ 

Differentiating (6.11) with respect to J, a straightforward calculation yields 

M(0,1)W^„ = 0, (614) 
q J 

for e — > 0, J — > and A 9 /i = 1, where M(e,h) = dfj(e,h)/dh is the mag- 
netization. Note that (dq/dJ) is independent on /i. From (6.14), we obtain 
(dq/dJ) J=0 = 0, and as a result 

= 0. (6.15) 



§ 7 The nonvanishing derivatives 

In this section, we show the non-universal behaviour of the relevant system 
when it satisfies the conditions 1) and 2). We use the formulas (5.3) and 
(5.6) derived from 1) and 2), respectively, and (6.10) and (6.15) result from 
the symmetry of the model. The properties (5.6), (6.10) and (6.15) are valid 
when or (or both and n^) are odd. Let us consider the weighted 
susceptibility 

X = P»b S 9ioio ( S io s jo ), ( 7 - 1 ) 

iojo 
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where gi j is the sign corresponding to the emerging order as gi j = sgn(gj g JO ) 
All we have to do here is to differentiate (7.1) in terms of Jm and to show 
that the second dominant term shows the logarithmic singularity Xo log e with 
e = (T-T c )/T c . 

Let us introduce the following notations. The interaction Ti, ki is expressed 

as 

H kl = -Jj20 t , and O^Ofof, (7.2) 

i 

where Jm is written as J for simplicity. The vector denotes the coordinate 
of the spin Sj and Ri = {r^ , • • • , r in } denotes the coordinate of the spin prod- 
uct Oi — Sj 1 • • • s in . The expectation value of Oi is expressed as a function 
of Ri as 

(Oi) = c (R i )+c 1 (R i )e\oge. (7.3) 

We also write Ri C lZ m (or % G 7^ m ) when r s G 7£ m for all r s E Ri. 

Here we derive finite derivatives of the exponent 7. The first-order deriva- 
tive of (7.1) is 

(|j),=o = FA E 0** EK^oWh - (si s j0 ) (Ol k) Of\], (7.4) 

where ( ) denotes the expectation value taken for J = 0. All the terms 
cancel except the cases: 

1-1) Sj , Sj G 7^^ (or Sj , Sj G 7^;), and both and are even, and 

1-2) Si Q G 7ifc, Sj Q G 7i;, and both and are odd. 

The derivatives (d'j/dJ)^, (dT c /dJ)o of the latter case 1-2) vanish as shown 
in section 6. As a result, we have only to treat the case 1-1) as the first- 
order derivative. This case is a simple generalization of the argument for the 
eight-vertex model in section 4. The derivative (7.4) is written as 

(%j)j=o = «EfeE[(^o^of)o - (s i0 s j0 )(Ol k) ) }(O%. (7.5) 

k)jt> i 

This is a generalization of (4.3). From (3.2) and (5.3), we obtain the following 
temperature dependence 

(%j =0 ~ XoE^M^)- + a 1 (R i )\oge\[c (R i ) + Cl (R t )e\oge], (7.6) 
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where ao(Ri), ai(Ri), co(Ri) and c\(Ri) are some constants. The term Gi = 
Gi^Gf* is factorized so that the inside of the brackets [•••][•••] in (7.6) is 
positive. We generally cannot omit the term ai(Ri). The sum J2i Gi^ a o(Ri) 
and EiG^a^Ri) appear as the coefficient of {dxk/dJk) j k=0 where Xk is 
the susceptibility of the model described by the Hamiltonian H k = 7i k — 
Jk J2i Oi and they are finite in the thermodynamic limit. The coefficient of 
loge, namely 

E&M^M^) + a 1 (R i )c (R i )), (7.7) 

i 

is finite because ^ Gi^ ao(Ri) , J2i Gi^ cii(Ri) , Gi^co(Ri) and GPci(Ri) are all 
finite. This is the first order derivative of 7. 

For all the cases except when both and have even symmetry 
(i.e. the case 1-1: both and are even), we obtained from (6.10) and 
(6.15) 

( aj)o = and (^)o = 0. (7.8) 

In these cases, we have to show that the second-order derivatives are non- 
vanishing and finite. From (3.1) and (7.8), the second-order derivative of the 
susceptibility x is 

- W-7(0)^(|^)o - (f^)ologe(O)]. (7.9) 

Hence the existence of the logarithmic singularity is the sign of continuously 
varying critical exponents. All we have to do is to find a term proportional 
to x l°g e i n the second-order derivative of x and to show the coefficient of 
Xo log e is finite. 

Differentiating x — P^b J2( s i s jo) with respect to J twice, we obtain 

«oio i j 

(7.10) 

where we have used that or (or both and n®) are odd and the 
expectation value (s^ ■ • • s in ) equals zero when n is odd. The following case 
remains nonvanishing: 

II-l) s io , Sj E Hk (or s io , sj £ Hi), is even and is odd, and 
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II-2) Sj , Sj G Hk (or Sj , Sj G 7Yz), and both n( fe ) and are odd. 

Both cases can be treated simultaneously. From (3.2) and (5.3), the deriva- 
tive (7.10) shows the following temperature dependence 

Xo E E GiGj i a o(Ri, Rj)- + «i (-Ri, Rj) log e] [c (-Rj, + ci Rj)e log e] , 
» j e 

= XoT.T.&GM^Rj), (7.H) 

where -Rj) is positive and ao(Ri, Rj), ai{Ri, Rj), co(Ri, Rj) and c\{Ri, Rj) 
are some constants. This summation can be regrouped by {7Z m } and using 
the condition (5.6) (i.e. the condition 2: Y^ien m & = 0)> we obtain 

E E X! GiGj^(R m + AR mi , i? m / + AR m fj) 

mm' i€TZ m j€TZ m i 

- XoE E GiGjiARmi- \7)(AR rn , r VMR m ,R m ,), (7.12) 

mm' i&l m 3&^ m i 

where i? m C 7?. m and i? m / C 1Z m i are fixed for each 1Z m and 7£ m / , respectively, 
and we have used the following notation 

AR mi ■ V = Ar mih ■ V h H h Ar mnin • V in 

Ai^-V = Ar Sjl • v;, + • • • + Ar mUn ■ V jn (7.13) 

where Vj and are the gradient operating to the coordinate of Sj and A mi - V 
and A m ij • V' operate on the first and the second arguments of Q(R m , R m >)- 
The first summation in (7.12) is classified by the distance between 1Z m 
and 7Z m r as 

E=E E , (7-14) 

mm' r mm'\m—m'\=r 

where \m — m'\ = r denotes that r < min{|rj — 7>| | G lZ m ,ri' G 7^ m '} < 
r + Ar and Ar is a constatnt comparable to the mean size of {7Z m }- 

As is a smooth and decreasing function of r and all {1Z m } denote finite 
regions containing a finite number of spins, each term in (7.12) is bounded 
by 

QiQj(AR mi ■ V)(AR m ,j ■ V)Q < const. (AR ■ V)(AR ■ V)Q, (7.15) 
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where 

AR- V = Ar • + ••• + Ar • V in , (7.16) 

and a similar equation is defined for Ai? • V. 

Then the coefficient of Xo^-Oge in (7.12) is bounded as 

Y^i^Rmi ■ V)(Ai? m , j • V)F(r, {R,}) 

r 

/CO 
d 2 r(Ai?- V)(A J R-V')F(r, {R m }), 

(7.17) 

where 

F{r, {R l })= E QiQjMRuRj^iiRuR^+aiiR^R^coiR^Rj)} 

mm'\m—m'\=r i£lZ m 3 £Tt m i 

(7.18) 

is a finite function of r (see (7.7)). 

Our purpose is to show (7.17) is finite. The first sum is of course finite 
and the second term is bounded by 



r°° f Or Or 

A / d 2 rV 2 F(r, {R,}) = A / dS— ~ A • 27rr — | r=ro < oo, (7.19) 

Jro Jr=ro Or OT 

where A is some constant. 



§ 8 Generalization 

The condition 2) can be generalized in the following two points. 

At first, it is straightforwald to generalize the form of the interaction TCu 

as 

W« = ^d W« = -J^^O^OS, (8.1) 

where each 7ij^ satisfies the condition 2). 

Next, the values of and have been fixed. This condition is nec- 
essary for showing the cancellations of correlations in the zeroth- and the 
first-order terms in (7.12) using the condition J2ieK m Gi = 0. However, our 
argument is valid for more complicated interactions i.e. the case that 
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and depend on the region lZ m . Each contribution from each lZ m to the 
derivatives is classified into the cases shown in sections 6 and 7. As for the 
first-order derivatives, there is no difference. For the second-order derivative, 
all the contributions to (<9 2 7 /dJ 2 ) coming from lZ m and TZ m i vanish except 
when both n$ + n^j and n$ + n®, are even. (Otherwise the argument in 
section 6 is valid for each term and the contributions to the derivative is van- 
ishing.) The nonvanishing case can be treated in the same way as in section 
7. 

§ 9 Conclusion 

We have derived that some Ising systems satisfying the conditions 1) and 2) 
in section 2 show non-universal critical behaviour. The perturbational expan- 
sion in terms of the interaction J is performed. This method results in the 
exact first-order derivatives (4.6) in the case of the eight- vertex model. We 
have used the conditions (5.3) and (5.6), which is derived from the conditions 
1) and 2), respectively, and (6.10) and (6.15) which result from the symme- 
try of the model. The existence of finite and non-zero derivatives {d^/dJ)® 
or ((9 2 7/<9J 2 )o is an evidence of continuous variation of critical exponents. 
These derivatives are derived in (7.4)-(7.7) and (7.10)-(7.19) for the first- 
and the second-order derivatives, respectively. Finally, the straightforward 
generalization of the condition 2) is commented in section 8. This condition 
is valid for generalized spin S Ising models or can easily be generalized for 
other classical systems. 
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